Journal of Research'of the National Bureau of Standards 



Vol. 58, No. 6, June 1957 



Research Paper 2764 



Thermal Design of Large Storage Vessels for 
| Liquid Hydrogen and Helium 1 

Russell B. Scott 

Principles of insulating vessels for the storage and transportation of liqueried gases are 
discussed. The use of shields cooled by the vapor leaving the liquid being stored is analyzed, 
and formulas, graphs, and tables are presented for the convenience of the designer of large 
vessels for storing and transporting liquid hydrogen and liquid helium. An analysis is 
given of the use of escaping vapor to intercept heat entering by conduction through solid 
supporting members. The evidence presented indicates that substantial improvements are 
possible in the design of such vessels. 



1. Introduction 

The current expansion in large-scale production 
of liquid hydrogen and helium has increased the 
importance of efficient storage and transport vessels 
for these expensive liquids. Published information 
on the design of low-loss containers includes descrip- 
tions [1,2,3,4,5,6] 2 of small laboratory vessels de- 
signed for storage of liquid helium. These, of course, 
are suitable also for storing liquid hydrogen. Some 
large, rugged, transportable liquid-hydrogen vessels 
have been constructed but have not yet been de- 
scribed in the open literature. It is the purpose of 
this paper to discuss some of the general principles 
relating to the thermal design of low-loss vessels and 
to present analyses of methods that can be used to 
improve the performance of an insulated vessel by 
utilizing the refrigerative effect of the escaping vapor. 
The material presented will relate specifically to the 
design of large vessels for liquid hydrogen and helium, 
although the principles used can be employed in the 
design of vessels for other liquefied gases. 

The invention of the vacuum-insulated vessel for 
liquefied gases by James Dewar in 1892 was a break- 
through in the field of thermal insulation that has 
not yet been matched by further developments. All 
of the advances since Dewar's time have been im- 
provements on Dewar's original concept, usually 
means of reducing radiant-heat transfer by attaining 
surfaces of higher reflectivity or by interposing shields 
that reflect or intercept the radiant energy. The 
Dewar vessel is illustrated in figure 1. At its best, 
the vacuum is so good that heat transfer by residual 
gas is almost negligible. Accordingly, the designer 
is mainly concerned with reducing the heat trans- 
ferred by mechanical supports or conveyed across the 
insulating vacuum by radiation. The heat trans- 
ferred by the supports is a problem involving the 
mechanical strength required for a specific applica- 
tion and has no general solution. If there is no 
restriction on the dimensions of the vessel, the heat 
conducted through the supports can usually be made 
quite small simply by using supports of low thermal 
conductivity and making their lengths great enough. 

1 This work was supported by the U. S. Atomic Energy Commission. 

2 Figures in brackets indicate* the literature references at the end of this paper. 



Even in restricted space the clever designer can 
usually find a way to provide a long thermal path 
through the solid supports. On the other hand, the 
heat transferred by radiation is little affected by the 
thickness of the insulating space. However, this 
latter mode of heat transfer can be reduced by inter- 
posing thermal shields between the warm and cold 
surfaces of the vessel. The shield is an opaque, 
highly reflecting surface, geometrically similar to the 
surface of the liquid container and suspended with 
the minimum of thermal contact approximately mid- 
way between the inner and outer walls of the vacuum 
space. If the shield has the same reflectivity as that 
of the boundary walls and contains no openings or 
irregularities, the radiant-heal transfer will be re- 



FILL AND VENT TUBE, THIN, 
POORLY CONDUCTING WALL 



HIGH VACUUM 




CHARCOAL 
ABSORBENT 
CHA 



STORED LIQUID 



Figuke 1. Metal Dewar vessel. 

Commercial vessels of this type usually are mounted in an outer protective 
shell of sheet steel. 



317 



duced by about one-half. Additional shields will 
further reduce the heat transferred by radiation. 
However, the design of mechanical supports, which 
will hold the shields in place without providing objec- 
tionable heat leaks, becomes a formidable task as 
the number of shields is increased. Also, openings 
and irregularities that are needed to accommodate 
the piping and supports for the inner container 
greatly reduce the effectiveness of this type of 
multiple shielding. 

2. Evacuated-Powder Insulation 

Fortunately, a method has been discovered that 
does not complicate the mechanical design and yet 
provides a barrier to thermal radiation comparable 
with multiple shields. This method consists of filling 
the insulating vacuum space with a fine light powder, 
such as perlite, siLca-aerogel, carbon black, or di- 
atomaceous earth. If the thickness of the insulating 
space is sufficient, filling with powder will greatly 
reduce the radiant-heat transfer. Moreover, the 
vacuum requirement is quite modest; the minimum 
heat transfer is very closely approached when the 
interstitial gas pressure is reduced to 10 \x of mercury 
in an insulating powder (e.g., perlite or silica-aerogel) 
with one boundary at room temperature, 300° K, 
and the other at the temperature of liquid nitrogen, 
77° K. A layer of evacuated powder several centi- 
meters thick will effectively stop radiation and yet 
not provide an objectionable amount of thermal 
conduction through the touching particles of powder. 
When the space between the walls must be small, of 
the order of 1 cm, the thermal conduction through 
the powder is objectionable, so in this case the most 
effective insulation is a high vacuum bounded by 
highly reflecting walls. Because the principal func- 
tion of the powder is to stop radiation, it is useful 
only when radiation constitutes an important heat 
leak. If an insulating vacuum is bounded by highly 
reflecting surfaces, such as clean annealed copper, 
silver, or aluminum, and the surfaces are at 77° and 
20° K, respectively, adding powder may actually in- 
crease the heat transfer by providing paths of solid 
conduction. The radiant-heat transfer between the 
surfaces at these temperatures is so small that its 
elimination by adding powder is of little conse- 
quence. Radiation from a surface at 77° K is only 
about 0.004 of that from an identical surface at 300° 
K, because it is proportional to the fourth power of 
the absolute temperature. Of course the thickness 
of the insulating space must be considered in such a 
comparison. It is probably more correct to say 
that for boundaries at 77° and 20° K a rather large 
thickness of powder would be required to compete 
with high-vacuum insulation. 

It appears then that evacuated-powder insulation 
is very desirable if the warm boundary is at ordi- 
nary room temperature and the insulating space is 
thick enough (10 cm or more) so that solid conduc- 
tion through the powder is not serious. 

Of course, evacuated-powder insulation is not new. 
It has been used for many years in large industrial 
containers for liquid oxygen and liquid nitrogen. It 



is particularly suitable for large containers. Now 
that liquid hydrogen is being stored and transported 
in large quantities, evacuated powder appears at- 
tractive as the insulation to use in large liquid- 
hydrogen vessels. Evacuated powder has been 
used in conjunction with a liquid-nitrogen-cooled 
shield for some vessels holding 750 liters of liquid 
hydrogen, and there are now under construction some 
liquid-hydrogen vessels of 6,000-liter capacity that 
utilize only powder insulation (no refrigerated shield) . 

3. Dependence of the Rate of Evaporation 
Upon the Size and Shape of the Vessel 

It is quite easy to formulate for liquefied-gas con- 
tainers some general relationships between the size 
and shape of the vessel and the rate of evaporation of 
the liquid it contains. First, because the heat leak 
into the container increases as the surface area 
increases, the most favorable shape is the sphere 
(the greatest ratio of volume to surface). The 
sphere also has some advantages of mechanical 
strength but in large sizes it may be expensive to 
fabricate. The next most favorable common shape 
is the square cylinder in which the length equals the 
diameter. The practical approach to this shape is 
the cylinder with standard commercial dished or 
elliptical heads. This design has good mechanical 
strength and is readily fabricated. It is stronger 
and has a somewhat larger ratio of volume to surface 
than the cylinder with the flat ends. However, in 
many cases, particularly vessels that are to be trans- 
ported by rail or highway, the maximum diameter is 
fixed. 

Thus most of the designs for large insulated tanks 
for storing and transporting liquefied gases can be 
assigned to two categories: (1) fixed shape and (2) 
fixed diameter. It is of interest to see how the per- 
formance of such containers varies with the size of 
the vessel. The effectiveness of a container for 
liquefied gas is commonly given as the fraction (or 
percent) of the total contents that evaporates in 
1 day. 

3.1. Vacuum-Insulated Vessels 

For vacuum-insulated containers of fixed shape, 
that is, containers of geometric similarity but of dif- 
ferent sizes, the heat leak by radiation and residual- 
gas conduction is proportional to the surface area of 
the inner container, that is, proportional to L 2 , where 
L is a linear dimension of the container. The heat 
leak through the supports is also proportional to Z 2 , 
because the supports will have a length proportional 
to the linear dimension and a total cross-sectional 
area proportional to the volume (and weight) of 
the inner container. Then the total rate of heat leak 
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where C\ is a constant and V is the volume of the 
inner container. The fractional rate of evaporation 

is 

w _ (7 1 v r(2/8) _ a 
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where L v is the heat of vaporization of the liquid 
and p its density. 

3.2. Vessels With Porous Insulation 

For containers insulated with porous materials, 
e. g., evacuated powders, the situation is somewhat 
different because, for similar geometry, the thick- 
ness, and therefore the effectiveness of the insulation, 
is proportional to the linear dimension of the con- 
tainer. The relation between the size and the 
performance of such vessels is not simple unless the 
heat conducted by the supports can be neglected. 
In this case the total heat leak is proportional to 
the linear dimension, so the fractional evaporation 
is proportional to l/V (2/3) . Thus it is seen that for 
vessels utilizing this type of insulation, large size is 
even more advantageous than it is for high-vacuum- 
insulated vessels. 

3.3. Vessels of Fixed Diameter and Different Lengths 

The case of the container for rail or road trans- 
portation, in which the maximum diameter cannot 
exceed a certain value but the length may be varied, 
can be conveniently analyzed on the assumption that 
the ends are spherical, so its shortest length and 
minimum volume are those of a sphere having the 
permitted diameter. In this case the thickness 
of the insulation is greatly restricted, so it will be 
assumed that the thickness is constant. This 
means that the heat leak through the insulation 
will be proportional to the surface area both for 
high-vacuum insulation and for porous insulation. 
It ^should be emphasized that this does not imply 
that the effectiveness of the various insulations is the 
same. It only means that in this circumstance 
the heat leak for unit area is constant for a given 
type of insulation and is not a function of the total 
area. Then, if we again ignore the heat leak through 
the supports, 

W=C [47rr 2 +27rr(i-2r)], 

where C is a constant, r is the maximum permissible 
radius, and L the length of the container. The 
fractional evaporation is 



W 



2CL 



L vP V L vP r(L-*r) 



As one would expect, the fractional evaporation 
changes very little with length, decreasing by only 
one-third when the length of the container is changed 
from 2r to infinity. 

3.4. Comparison of Evaporation Rate Versus Size for 
Different Vessels 

Figure 2 shows graphically how the fractional 
rate of evaporation depends upon the volume of 
the container for the three conditions just discussed. 
It should be pointed out that the three curves should 
be compared with each other in regard to trend, 
not in regard to magnitude. The actual rates of 
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Figure 2. Dependence of the evaporation rate upon the size 
of the vessel. 

A, Vessels of fixed diameter and various lengths; B, high-vacuum-insulated 
containers or containers having constant insulating thickness and similar geom 
etry; C, powder-insulated vessels of similar geometry throughout, including a 
constant ratio of insulating thickness to the other linear dimensions of the vessel. 

evaporation can be determined only by computing 
the heat leaks for the particular insulation being 
used. 

3.5. Convenient Formula for Porous Insulation 

A good approximation for the heat transfer 
through a porous insulation of constant thickness 
in a vessel of conservative shape is given by 

where k is the mean effective thermal conductivity 
between the tempertaure T 2 and T x , T 2 is the tem- 
perature of the outer surface, T\ is the temperature 
of the inner surface, A x and A 2 are the areas of the 
inner and outer surfaces respectively, and t is the 
thickness of the insulation. This formula is exact 
for concentric spheres. For cylinders with spherical, 
dished, or elliptical heads it has been estimated that 
the formula will give results within 5 percent if the 
thickness of the insulating space is not more than 50 
percent of the radius of the cylinder. It is cor- 
respondingly reliable for other shapes in which the 
thickness of the insulation is a small fraction of the 
minor perimeter of the vessel. In reporting values 
of the thermal properties of insulators, the mean 
effective thermal conductivity is often given between 
commonly encountered temperatures. If the avail- 
able data are presented in the more general form 
of instantaneous conductivities as a function of 
temperature, the mean effective thermal conductivity 
can be obtained from the relation 
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where k is the thermal conductivity at temperature 
T. If k is a linear function of T, k is the arithmetic 
mean conductivity between T x and T 2 . 



421239—57- 



319 



THIN STAINLESS STEEL 



LIQUID NITROGEN 



CHARCOAL 
ADSORBENT 



Figure 3. 




r COPPER 



HIGH VACUUM 



LIQUID HELIUM 
OR HYDROGEN 



Simplified cross section of a Dewar vessel with a 
liquid-nitrogen-cooled shield. 



4. Shield Cooled With Liquid Nitrogen 

The vessels described in references [1] to [6] use 
liquid nitrogen to reduce the rate at which heat 
reaches the liquid helium or hydrogen being stored. 
In some cases the Dewar is simply immersed in liquid 
nitrogen contained in a larger vessel. Other designs 
are of unit construction, like that shown in figure 3, 
in which there is a liquid-nitrogen-cooled shield sus- 
pended in the vacuum between the inner container 
and the warm wall. Most of the heat entering from 
the exterior is absorbed by the relatively inexpensive 
liquid nitrogen. This design is the most effective 
yet produced and is the most practical for storing 
small quantities of liquid hydrogen and helium. 
However, as the size of the vessel increases, the frac- 
tional evaporation rate without the shield decreases, 
and eventually a point is reached at which the extra 
initial cost of the liquid-nitrogen-cooled shield and 
the cost and bother of supplying liquid nitrogen are 
not worth the extra saving of the liquid being stored. 
This change-over point depends on the relative costs 
and heats of vaporization of liquid nitrogen and the 
liquid being stored. 

5. Vapor-Refrigerated Shields 

The vapor leaving the inner container of an in- 
sulated vessel is normally at a temperature very near 
the- boiling temperature of the liquid being stored. 
This cold vapor constitutes a source of refrigeration 



that can be used to absorb some of the heat leaking 
into the vessel. In most vessels some of this refriger- 
ation is utilized more or less accidentally in absorb- 
ing some of the heat which flows by solid thermal 
conduction down the vent pipe. An analysis of this 
process will be given later. 

A way to make good use of the refrigeration avail- 
able in the escaping vapor is to have it cool a shield 
suspended between the inner liquid container and 
the outer warm wall of the vessel. This device can 
be applied to high-vacuum-insulated vessels, as well 
as vessels utilizing evacuated-powder insulation. The 
openings in the shield required to admit piping and 
supports or reentrant irregularities that may be pro- 
duced by attaching piping to one side of a shield will 
have little effect when powder insulation is used, 
but when the shield is in an empty space these open- 
ings and irregularities may somewhat reduce the 
effective over-all reflectivity of the shield. More- 
over, in a vessel having sufficient room between the 
walls to accommodate a shield, the insulation thick- 
ness may be great enough to show a substantial im- 
provement when powder is used. Because of the 
simpler construction and the fact that the ideal is 
more nearly approached, it is believed that the vapor- 
cooled shield is more aptly applied to a vessel with 
evacuated-powder insulation. 

Some commercial powder-insulated vessels make 
partial use of the cold escaping vapor by having it 
flow through several turns of pipe that lie between 
the inner and outer walls. It is believed, however, 
that the available refrigeration can be used more 
effectively if it is made to cool a shield located at 
the optimum position between the inner and outer 
walls of the insulating space. 

The following analysis presents the general method 
of determining the optimum location of a vapor- 
cooled shield and the expected gain in performance. 
The results are applied to the specific case of spherical 
containers for liquid hydrogen and helium with evac- 
uated perlite insulation. 

With evacuated-powder insulation filling the entire 
insulating space, the rate at which heat reaches the 
liquid is 



r 2 —r l 



(1) 



where k 1 is the mean thermal conductivity of the 
powder in the temperature range T x to T 2 , r± is the 
radius of the liquid container, r 2 is the radius of the 
vapor-cooled shield, m is the mass rate of evapora- 
tion, and L v is the heat of vaporization of unit mass 
of the liquid. 

The rate at which heat reaches the shield is 



W 2 = 



^k 2 r 2 r z {T z -T 2 ) 



n—r 2 



^mL^rhAh, 



(2) 



where r 3 is the radius of the outer boundary of the 
insulation, k 2 is the mean thermal conductivity of the 
powder in the temperature range T 2 to T 3 , and AA 
is the change of specific enthalpy of the vapor as it 
is warmed froM T Y to T 2 . 
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Dividing (1) by (2) gives 



L v = ki (r z —r 2 ) n (T 2 —T 1 ) i 
L v +Ah k 2 (r 2 — ri) r 3 (T z —T 2 ) 

Let r z =ari and r 2 =br h then 

q(&-i) _fa z g +Afe c^ -ro 



From (1), 



isr*5=i< r '- r '>- 



(3) 



Because 6 is a function of a and # only, and g and & 
are functions of T 2 only, m/ri is a function of a and 
T 2 only. 

For each value of a there is an optimum value of 
T 2 and a corresponding optimum value of b, which 
will make m\r x a minimum. Because, for a given 
insulation, k u k 2j and g depend only upon the proper- 
ties of the liquid and vapor, the fixed boundary 
temperatures T t and T z , and the variable tempera- 
ture T 2 , we can compute values of g for a series of 
temperatures T 2 in the range 2Y\ T 2 < C T z and use 
these g's for all subsequent computations dealing 
with the particular liquid being stored. 

Then, for a given value of a, values of rh/ri can 
be computed and the minimum can be determined 
by graphical means. By plotting m\r x versus the 
position of the shield, (b — l)/(a— 1), the location of 
the minimum will also yield the optimum value 
of b, the important design parameter. It will be 
noted also that for this vessel the rate of evaporation 
is proportional to the first power of the radius or to 
(volume) 173 , if heat conduction through supports is 
neglected. Thus this design offers the same premium 
for large size as that shown in curve C of figure 2. 

5.1. Liquid-Hydrogen Containers 

A series of computations was made for values of a 
ranging from 1.05 to 1.5 for a vessel containing 
liquid parahydrogen with evacuated-perlite insula- 
tion. The thermal conductivity of perlite is shown 
in figure 4, a curve based on two determinations 
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Figure 4. Thermal conductivity of — 80-mesh perlite. 



made at NBS, one giving the mean effective con- 
ductivity between : > ,()4° and 76° K, the other the 
mean conductivity between 70° and 20° K. These 
data arc admittedly not as complete as they should 
be for the purpose at hand, but they will serve bo 
demonstrate the usefulness of the vapor-cooled shield 
in this application. Errors arising from the linear 
interpolation used in figure 4 are not expected to 
exceed 20 percent. Additional measurements of the 
thermal conductivities of powders as functions of 
temperature are in progress at NBS. The values of 
Ah were taken from the paper of Wooley, et al. [7]. 

Figure 5 shows the results of a representative 
computation for a liquid-hydrogen container. It is 
seen that at the optimum position of the shield, 
about 35 percent of the distance from the cold wall 
to the warm wall, the shielded vessel will have an 
evaporation rate of about 38 percent of the un- 
shielded vessel. It is seen that the optimum location 
of the shield is not critical; the curve of m/m* has a 
very flat minimum. It was found that the effective- 
ness of the shield, ra/m*, at its optimum location is 
practically the same for values of a ranging from 
1.05 to 1.5. 

It will be noted that the shield is somewhat 
warmer than liquid nitrogen. The rate of evapora- 
tion will be lower if a liquid-nitrogen-cooled shield is 
used. The decision about the use of this method 
rather than liquid-nitrogen shielding will depend 
upon the relative costs of liquid nitrogen and liquid 
hydrogen and the cost of keeping the liquid nitrogen 
shield replenished. Also, the complete design study 
will include heat leaks from piping and supports, 
which cannot be resolved here because they are 
peculiar to the specific application and mechanical 
design. 

It is remarkable also that the temperature of the 
shield and its optimum location, that is, its distance 
from the inner container divided by the distance 
between the inner and outer walls, (6 — l)/(a— 1), are 
practically independent of the value of a. This is a 
consequence of the temperature dependence of the 
thermal conductivity of the powder and the specific 
heat of the vapor. As more complete data on the 
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Figure 5. Optimum position and temperature of a vapor- 
cooled shield in a powder-insulated vessel for liquid hydrogen. 
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thermal conductivity of powder insulations become 
available, the optimum position should be recom- 
puted. 

To gain an idea of the performance that might be 
realized in a practical vessel, the results of these 
computations were applied to a spherical container 
for liquid hydrogen having an inner shell 2 m in 
outside diameter and an outer shell 3 m in inside 
diameter. This will have a capacity of slightly more 
than 4,000 liters. The evaporation rate resulting 
from heat leak through evacuated-perlite insulation 
was found to be sufficient to evaporate 25 liters, or 
0.62 percent, of the contents per day. A vapor- 
cooled shield reduces this to 0.24 percent per day. 

5.2. Containers for Liquid Helium 

It is quite apparent that the effectiveness of the 
vapor-cooled shield depends upon the ratio of the 
heat of vaporization of the liquid being stored to the 
specific heat of the vapor. Accordingly the appli- 
cation of this device to the storage of liquid helium 
is very attractive because of the favorable ratio 
between the specific heat of the vapor and the heat 
of vaporization. 

The results of computations for a =1.5 fully 
justified this optimism because it was learned that 
the vapor-cooled shield of a powder-insulated con- 
tainer for liquid helium at its optimum location will 
have a temperature of 62° K, substantially lower 
than that of a liquid-nitrogen-cooled shield. The 
optimum position of such a shield was found to be 
(b— I) /(a— 1)=0.25. The designer should again be 
cautioned against using this value after more data 
on the thermal conductivity become available. 
Measurements now in progress will undoubtedly 
yield more accurate data for this application. 

Applying these results to the 4,000-liter spherical 
container described above, it was found that the heat 
leak through the insulation would cause liquid helium 
to evaporate at a rate of 7.7 percent per day from 
the unshielded vessel and 0.7 percent per day from 
the shielded vessel. 

5.3. Combination of High- Vacuum and Evacuated- 
Powder Insulation 

A logical refinement of the foregoing design is to use 
a high vacuum in the space between the coldest wall 
and the vapor-cooled shield, and powder insulation 
in the outer space. If the surfaces bounding the 
vacuum space have a low emissivity, this arrange- 
ment will result in a considerable improvement, 
gained at the expense of greater construction costs. 
The heat leak to the inner vessel, exclusive of solid 
conduction, can be computed by a method rather 
similar to that used earlier. In this case, however, 
the thickness of the high-vacuum space should be 
made as small as practical, so we can idealize the 
problem by assuming that the vacuum space has 
negligible thickness (r 1 ^r 2 ) . Then the rate at which 
heat is transmitted to the inner container by radia- 



tion from a shield at the temperature T 2 



is 



2-e 



±Tii(Ti-Tl)=mL„ 



(4) 



where a is the Stefan-Boltzman radiation constant, 
e is the emissivity of the surfaces facing the vacuum 
space, ri is the radius of the inner container, T 2 and 
T\ are the temperatures of the shield and inner con- 
tainer respectively, m is the mass rate of evaporation 
of the liquid being stored, and L v is its heat of vapor- 
ization. The heat that reaches the shield by con- 
duction through the powder is 



kAirrMTt-Ti) 



(rz—n) 



=mL v -\-mAh, 



(5) 



where h 2 is the mean thermal conductivity of the 
powder, r 3 is the radius of the outer shell at tempera- 
ture T 3 , and Ah is the change in enthalpy of the 
vapor as it is warmed from T\ to T 2 . Dividing eq 
(4) by (5) and letting a=r 3 /ri gives 



aenia-lXTI-Tt) 



L v +Ah (2-e)k 2 a(T z ~T 2 ) 



(6) 



Because Ah and k 2 are functions of T 2 that are not 
available from analytic expressions, the simplest 
method of solving this equation for T 2 is to use suc- 
cessive approximations. A general solution cannot 
be given. The equation must be solved for specific 
values of r x and a. For this reason the performance 
of the 4,000-liter vessel with an inner diameter of 2 
m and an outer diameter of 3 m, described earlier, 
was computed so that the results can be compared 
with those resulting from other insulation methods. 
It was found that with liquid hydrogen, the shield 
temperature will be 129° K and the evaporation rate 
will be 0.13 percent per day. When this system is 
used with liquid helium, the temperature of the 
shield is 86° K and the evaporation rate 0.33 percent 
per day. These computations assume an emissivity 
of 0.02 for the surfaces facing the vacuum space and 
ignore the heat leak through the supporting members. 

5.4. Discussion 

Table 1 shows the rates of heat transfer in terms of 
loss of liquid hydrogen and liquid helium in percent 
per day for several different types of insulation for a 
vessel of 4,000-liter capacity. These results do not 
take account of heat transfer by conduction through 
the supports nor heat conduction through the resid- 
ual gas in a high vacuum. It should be emphasized 
that the designer and fabricator have a great amount 
of control over these latter sources of heat leak. 
The length of thermal path through insulating sup- 
ports is particularly pliant to the ingenuity of the 
designer. Figure 6 gives a suggestion for increasing 
the length of supports. 

The computations show that the effectiveness of 
the vapor-cooled shield in powder insulation, ra/m*, 
is practically independent of the thickness of the 
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Table 1. — A comparison of the heat leaks through various insulations for a 4,000-liter spherical Dewar 

The outside diameter of the inner vessel is 2 m and, when powder insulation is used, the outer shell has an inside diameter of 3 m. For the cases in which no 
powder is used, the diameter of the outer shell will be smaller. The emissivity of the surfaces facing the high-vacuum spaces is taken as 0.02. 



Insulation 


Evaporation rate caused by the heat flow through the insulation 
(excluding conduction by solid supports and assuming no con- 
duction by residual gas) 


Liquid hydrogen 


Liquid helium 


High vacuum a ._. - -- . - .._.._. 


%/Day 
4.0 

0.62 

.24 
(shield temperature =124° K) 

.13 
(shield temperature =129° K) 

.017 


%/Day 
49.0 

7.7 

0.70 
(shield temperature = 02° K) 

.38 
(shield temperature = 86° K) 

.21 


Evacuated perlite ._. _. ___ 


Evacuated perlite plus vapor-cooled shield ._ 


Evacuated perlite plus high vacuum a plus vapor-cooled shield __. _ _ 


Uigh vacuum a plus liquid N2 shield (77° K) 





SUPPORT ROD 
OR CABLE 



INSULATING 
SPACE 




RE-ENTRANT 
TUBE 



Figure 6. 



Reduction of thermal conduction through, supports 
by increasing their length. 



» These are idealized conditions. The vacuum space is assumed to be negligibly thin and the vacuum is assumed to be so good that heat conducted by residual 
gas can be ignored. 

thoroughly condensed when the container holds 
liquid helium. On the other hand, hydrogen is 
dissolved in, and is slowly released by, nearly all 
materials used to construct containers for cryogenic 
liquids, so provision should be made for its disposal 
if a high vacuum is to be maintained in vessels for 
other than helium. 

A very obvious extension is the use of multiple 
shields. Because the cold vapor has a substantial 
rcfrigerative value after it leaves the shield discussed 
here, it is apparent that it could be utilized to cool 
additional shields. 

Another source of refrigeration that could be 
utilized to cool the shield of a container for liquid 
parahydrogen is the heat absorbed by para-to-ortho 
conversion. If the proper catalyst is put into the 
vent tubes that cool the chield, the endothermic 
conversion to a higher ortho concentration will pro- 
vide extra refrigeration. The cooling produced by 
converting hydrogen of 99.8 percent para (the equi- 
librium concentration at the normal boiling temper- 
ature) to 38.6 percent para (the equilibrium concen- 
tration at 100° K) is 44 percent of the cooling pro- 
duced by warming the gas from the hydrogen boiling 
point to 100° K. 

Finally, a word should be said about the storage 
of liquid oxygen and nitrogen. Because the ratio 
of the heat of vaporization to the change in enthalpy 
of the vapor upon being warmed to room temper- 
ature is much greater for these materials than it is 
for hydrogen and helium, it should be expected that 
the value of a vapor-cooled shield will be correspond- 
ingly less. However, the principles just discussed 
are equally valid for this application and should be 
given adequate consideration in the design of vessels 
for the long-term storage of liquid oxygen and 
nitrogen. 



insulation, here measured by the constant a. The 
optimum location of the shield is slightly affected by 
the thickness of the powder. For the liquid-hydro- 
gen vessel, the optimum position varied from 39 per- 
cent of the way from the cold to the warm surface for 
a =1.05 to 32 percent of the way for a =1.5. The 
two computations for helium gave 30 percent for 
a=\.2 and 25 percent for a=1.5. These variations 
are rather trivial considering the flatness of the 
minimum in figure 5. 

The attainment of a first-rate high vacuum in an 
insulating space protecting the inner container 
requires first of all a completely leak-free construc- 
tion, and also a thorough out gassing of all the metal 
surfaces bounding the high vacuum. The use of an 
adsorbent cooled to the temperature of the liquid 
being stored is highly recommended. The attain- 
ment of a superior vacuum is more easily accom- 
plished in vessels for liquid helium storage than in 
liquid-hydrogen vessels because liquid helium pro- 
vides a temperature at which all gases except helium 
have a negligible vapor pressure. Helium is not 
likely to be found in an insulating vacuum unless 
there is a leak. It has been observed that com- 
mercial containers with liq 11 id-nitrogen-cooled shields 
have a much smaller heat leak when they are used 
to store liquid helium than when used with liquid 
hydrogen, probably because hydrogen gas is very 



6. Vapor-Cooled Vent Tube 

In those cases in which heat conduction through 
the supports or fill and vent tubes is a substantial 
part of the total, it may be worthwhile to use the 
escaping vapor to intercept part of the heat. This 
is an important consideration in the vent tubes of 
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some Dewars. The escaping vapor maintains good 
thermal contact with the wall of the tube and ab- 
sorbs a great deal of the heat that would otherwise 
enter and cause evaporation of the liquid being 
stored. The maximum possible saving from this 
process can be computed upon the assumption that 
the heat transfer between the vapor and the tube is 
perfect; that both have the same temperature at 
each level. Also, it is assumed that there is no 
lateral conduction or radiation to or from the tube. 
Because, in the temperature regions that are of in- 
terest, the thermal conductivity of the tube or sup- 
port is a function of temperature, the practical solu- 
tion of the problem can be simplified by assuming 
that this function is linear, that 

k=k +a(T-To), 

where k is the thermal conductivity at the variable 
temperature T, k is the thermal conductivity at the 
cold end, which lias the temperature T Qj and a is a 
constant. Then 



W=A[k +a(T- T )](dT/dx), 



(7) 



which is the equation for heat conduction where W 
is the heat current at any point in the tube, A is the 
cross-sectional area of the tube, and x is the dis- 
tance from the cold end. Now the issuing vapor 
will absorb some of the heat traveling down the tube, 
so that 

W=W +mC p (T-T ), (8) 

where W is the heat current that reaches the cold 
end of the tube, m is the mass rate of flow of the 
issuing vapor and Cp is its specific heat. 

Upon combining eq (7) and (8) and integrating, 
it is found that 



L/A=[lf(mC P )] 



in 



Wo+mC p AT 
W 



) 



(9) 



where L is the length of the tube and AT the tem- 
perature difference between the warm and cold ends. 
The thermal conductivities of alloys commonly 
used for vent tubes of metal Dewars are rather com- 
plex functions of temperatures. However, eq (9) 
can be employed by dividing the total temperature 
interval into smaller intervals in which the thermal 
conductivity can be considered to be a linear func- 
tion of temperature. Then eq (9) is applied to the 
lowest temperature interval, ATi=T 2 —Ti, and the 
dependence of W upon L X \A for various values of m 
is determined; L x is the length of tube lying in the 
lowest temperature interval. Then values of heat 
current W' Q at the bottom of the next higher tem- 
perature interval are determined by the relation 

W'o^Wo+mCvATL 

These values of W' are then used in eq (9) to com- 



pute values of L 2 /A for the second temperature in- 
terval, and the process is repeated as required. 
Finally, total values of L/A are obtained by adding 
the individual values, L n jA. 

For a stainless-steel vent tube for hydrogen it is a 
sufficiently good approximation to use just two tem- 
perature intervals, from 20° to 80° K and from 
80° to 300° K. In figure 7 are plotted values of W 
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Figuke 7. Results of computations of heat current reaching 
the cold end of a solid conductor counter to a stream of cold 
vapor. 
A,m = l g/sec; B, 10" 1 g/sec; C, 10" 2 g/sec; D, 10-3 g / S ec; E, 10-4 g/sec, F, w=0. 
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as a function of L/A for stainless-steel vent tubes. 
It should be remembered that these are limiting 
values, which assume complete heat exchange be- 
tween the escaping vapor and the vent tube. In 
practical vessels, it may be necessary to provide 
baffles or other means for promoting heat transfer 
between the escaping vapor and the vent tube to 
approach the thermal equilibrium assumed for 
these computations. For hydrogen, which has a 
variable specific heat, mC p can be replaced by 
m(Ah/AT), where Ah is the change of specific en- 
thalpy for the temperature interval AT. Values of 
thermal conductivities of stainless steel were taken 
from the data compiled by Powell and Blanpied [8]. 

7. Conclusions 

The foregoing analyses show quite clearly that 
important benefits can be realized by utilizing the 
refrigerative value of the vapor which escapes from 
an insulated vessel for liquid hydrogen or helium. 
The problems were idealized in order to make their 
solution simpler; thus in many practical designs t he- 
values reported in this paper may be substantially 
altered by the mechanical requirements. However, 
there is some definite experimental evidence that 
these computed values are not unrealistic. A 50- 
liter commercial container for liquid helium, similar 
to that illustrated in figure 3, protected with a 
liquid-nitrogen-cooled shield, showed an evapora- 
tion loss of 0.35 liter per day. If all this heat leak 
is attributed to thermal radiation from the shield 
at liquid-nitrogen temperature, the emissivities of 
the surfaces facing the vacuum space (assumed 
equal) are found to be 0.015, which is an excellent 
value for copper, better than the 0.02 assumed in 
table 1. This result implies that the heat conduction 
through residual gas and through the supporting 
neck of the Dewar was negligible. From the dimen- 
sions of the Dewar neck, 15.9 mm o.d., 0.5 mm 
thick, and 31.8 cm long (L/A=13Q cm -1 ), figure 7 
shows that if no vapor were escaping, the heat flow 
down the neck would be about 0.022 w, more than 
twice that required to cause the observed evapora- 
tion rate. However, with the observed rate of 
evaporation, 0.35 liter per day, 5.1 X10~ 4 g/sec, an 
interpolation 3 of the curves of figure 7 shows that 

3 The spacing of the lower parts of the curves of figure 7 indicates that logarith- 
mic interpolation should be satisfactory in the region under consideration. 



the heat that actually reaches the cold end of the 
neck is less than 10" 4 w, a result consistent with the 
observations. Thus, it appears that a performance 
very near the idea] has been achieved in a relatively 
small vessel with a shield cooled with liquid nitrogen. 
Accordingly it should be entirely practical to ap- 
proach the ideal in large vessels having either liquid- 
nitrogen-cooled shields or vapor-cooled shields. 

In many applications the neck of the 50-liter 
helium container described above is inconveniently 
small. It is possible to increase the neck diameter 
substantially, keeping the wall thickness and length 
the same, without seriously increasing the heat flow 
to the liquid helium. A value of L/A of 80 cm -1 
will give a heat current of very nearly 10 -4 w, only 
1 percent of the total heat leak. Doubling the 
diameter of the neck, making L/A=65 cm -1 , will 
result in a heat current at the bottom of the neck 
of 2X10 -3 w, a 20-percent increase. When the 
larger neck tubes are used, removable baffles should 
be provided, both to increase heat transfer between 
the escaping vapor and the tube, and to stop radia- 
tion down the tube. 



The assistance of James A. Brennan in performing 
the tedious computations is gratefully acknowledged. 
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